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Abstract. With more and more exoplanets being detected, it is paid closer attention to whether there are lives 
outside solar system. We try to obtain habitable zones and the probability distribution of terrestrial planets in 
habitable zones around host stars. Using Eggleton's code, we calculate the evolution of stars with masses less 
than 4.00 Mq. We also use the fitting formulae of stellar luminosity and radius, the boundary flux of habitable 
zones, the distribution of semimajor axis and mass of planets and the initial mass function of stars. We obtain the 
7^ ■ luminosity and radius of stars with masses from 0.08 to 4.00 Mq, and calculate the habitable zones of host stars, 

affected by stellar effective temperature. We achieve the probability distribution of terrestrial planets in habitable 
Q , zones around host stars. We also calculate that the number of terrestrial planets in habitable zones of host stars 

' is 45.5 billion, and the number of terrestrial planets in habitable zones around K type stars is the most, in the 

' Milky Way. 
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OO 1. Introduction Typically, stellar HZ is defined as region near the host 

star, where water at the surface of a terrestrial planet 

rn • Nowadays, more and more exoplanets have being indi- is in liquid phase (eg.. Hart, 1978; Kasting et al., 1993; 

rectly detected, and they were even directly detected last Franck et al., 2000; Noble et al., 2002; Jones et al., 2006). 

rn ; year (Marois et al, 2008). And theories about planetary Previously, it is thought that the boundary flux of HZ de- 

O ; evolution point out that most of stars own planets. It is pends on luminosity (eg., Kasting et al., 1993). However, 

^ ■ paid closer attention to whether there are lives around it ig pointed out that this flux also depends on stellar ef- 

other stars, not only by astronomy filed, but also by many festive temperature T^s (Forget & Pierrehumbert, 1997; 

other realms. A planet which is suited to life survive must wilhams & Kasting, 1997; Mischna et al., 2000; Jones 

satisfy two factors, one is that it is a terrestrial planet, and 2004; Jones et al. 2006). As the greenhouse effect can 

^ I the other is that it must be in the habitable zone (HZ) of ^aise the mean temperature of terrestrial planets (Clube 

. . , its host star. gt ^1., 1996), and the greenhouse effect is different for dif- 

Terrestrial planets are rocky planets from one to ten ferent radiation. For example, the lower Tcff, the more 

Earth masses with the same chemical and mineral compo- the infrared fraction in luminosity, and the more this frac- 

sition as the Earth (Valencia et al., 2006). However, others tion, the more the greenhouse effect for a given stellar flux 

argued that the mass of terrestrial planet could be down to (Jones et al., 2006). Thus, the distances at both inner and 

about 0.3 Earth mass, to retain its atmosphere over long outer HZ boundaries are farther to host star, with lower 

geological timescales and to sustain tectonic activity as Teff, than they would have been if the Tcs effect is not 

required for the carbon-silicate cycle to operate (Williams taken into consideration. And the distances are closer to 

et al, 1997; Menou & Tabachnik, 2003). host star with higher Toff. 
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It is well known that HZ widths for different stars are 
different. For example, the HZ width for a M type star is 
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only about one fifth to one fiftieth of the HZ width for a G 
type star (Tarter et al., 2007). However, it does not mean 
that the probability of terrestrial planet in HZ around a 
M type star is also about one fifth to one fiftieth of the 
probability around a G type star. Fortunately, planetary 
evolution has been calculated, which matched well with 
the observation data considering selective effects (Ida & 
Lin, 2004, 2005; Schlaufman et al., 2009). Ida & Lin (2005) 
gave the distribution of scmimajor axis and mass of plan- 
ets predicted by the Monte Carlo simulations. The masses 
of host stars are 0.20, 0.40, 0.60, 1.00 and 1.50 Mq, respec- 
tively. And this distribution can be used to calculated the 
probability of terrestrial planets in HZs around different 
host stars. 

Jones et al. (2006) gave the HZ flux at both inner and 
outer boundaries, as a function of T^s, but did not give 
the flux for host stars. Using Eggleton's code, we calcu- 
late stellar evolution and achieve the relationship between 
stellar mass and T^s- We get the HZ flux, at both in- 
ner and outer boundaries, for host stars at zero age main 
sequence (ZAMS) and at the terminal of main sequence 
(TMS). And the HZ flux for host star with mass 4.00 
Mq is about five times more than that of host star with 
mass 0.40 Mq at inner boundary and about four times 
at outer boundary. Then we obtain the HZs of host stars 
with masses from 0.08 to 4.00 Mq at ZAMS and at TMS, 
taken the Tgg effect into consideration. 

Using the distribution of semimajor axis and mass of 
planets (Ida & Lin, 2005), we obtain the probability dis- 
tribution of terrestrial planets in HZs around host stars. 
Then, we calculate that the number of terrestrial planets 
in habitable zones of host stars is 45.5 billion in the Milky 
Way. And the number of terrestrial planets in HZs for M, 
K, G and F type stars are respectively as 11.548 billion, 
12.930 billion, 7.622 billion and 5.556 billion, in the Milky 
Way. 

The outline of the paper is as follows: we present some 
input descriptions in Section 2, give our results in Section 
3, present some discussions in Section 4, and then finally 
in Section 5 we give our conclusions. 

2. Input descriptions 

2.1. The least lifetime of the main sequence for host 
stars 

Life emergence needs habitable environment with long 
enough time. Jones et al. (2006) thought that it needed at 
least 1.0 Gyr. And planets cover a presumed heavy bom- 
bardment phase as on Earth, at the first 700 Myr (Jones, 
2004; Lai, 2008). Hence, The lifetime of the main sequence 
(MS) for host stars is at least 1.7 Gyr, according to this 
issue. 

However, this requirement may be too strict. The iso- 
topic signature of carbon granules in the Isua formation 
(more than 3.7 Gyr) is interpreted as evidence for prokary- 
otic life (Rosing, 1999). And there is evidence that there 
was liquid water on Earth about 4.3 Gyr ago (Mojzsis et 



al., 2001). It is well known that the solar age is about 4.57 
Gyr (Bahcall et al., 1995). These mean that archaic life 
had come forth when our solar was about 800 Myr old, 
and the environment on Earth may be suited to biological 
evolution when solar was about 200-300 Myr old. Hence, 
we adopt that the least lifetime of MS for host stars is 
about 200 Myr. 

2.2. Input physics of stellar evolution 

We use the stellar evolution code of Eggleton (1971, 1972, 

1973), which has been updated with the latest input 
physics over the last three decades (Han et al., 1994; Pols 
et al., 1995, 1998). We set the convective overshooting pa- 
rameter, 5ov = 0.12 (Pols et al., 1997; Schroder et al., 
1997). We adopt metal mixture by Grevesse & Sauval 
(1998). 

We use OPAL high temperatures opacity tables 
(Iglesias & Rogers, 1996; Eldridge & Tout, 2004) in the 
range of 4.00 < log(T/K) < 8.70, and the new Wichita 
state low temperature molecular opacity tables (Ferguson 
et al., 2005) in the range of 3.00 < iog(r/K) < 4.00. 
And we have made the opacity tables match well with 
Eggleton's code (Chen & Tout, 2007; Guo et al., 2008). 

In our calculation, the value of metallicity is 0.02. As 
the lifetime of MS for star with mass 4.00 Mq is about 183 
Myr, we calculate the evolution of stars with masses less 
than 4.00 Mq. Stars with masses less than 0.50 Mq are 
spaced by AM = O.O5M0, and stars with masses more 
than 0.50 M© are spaced by AM = O.IOM©. 

2.3. The initial mass function of stars 

In the Milky Way, it is well known that the number of 
stars with different masses is different. In our calculation, 
we adopt the initial mass function of Kroupa et al. (1993). 

ciM-i-3, 0.08 < M < 0.5, 
C(M) = I C2M-2-2, 0.5 < M < 1.0, (1) 
C2M-2-7, 1.0<M<oo, 

where ^{M)dM is the probability of stars with masses 
between M and M-|-dM, in solar units. The mass function 
is normalized according to 

/ ^(M)dM = 1, (2) 

JO.08 

SO that ci = 0.255508 and C2 = 0.138704. 
3. Results 

3.1. Fitting formulae for stellar luminosity and radius 

The fitting formulae of stellar luminosity and radius have 
been completed by formers (Tout et al., 1996; Hurley et 
al., 2000). However, the studies before 2005 used the old 
Wichita state low temperature molecular opacity tables 
(Alexander & Ferguson, 1994), and there are some errors 
for the old opacity tables (Weiss et al., 2006). Therefore, 
it is required to fit these formulae over again. 
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Table 1. The fitting coefficients for Eqs. (3), (4), (5), (6) 
and (14), eg., 020 = 2.281557 x 10^. 
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3.1.1. Stellar luminosity and radius at ZAMS 

For stellar luminosity and radius at ZAMS, we fully use 
the former shapes (Tout et al., 1996), and achieve a series 
of new fitting coefficients. 



^ZAMS 



aiM^-s + a2Mi 



(3) 



.RzAMS = 

ai3 + ai4M2 + aisMS-S + M^S-S + aieMiS-^ ' 

where Lzams and i?zAMS are stellar luminosity and ra- 
dius at ZAMS, respectively, both in solar units. And M 
is stellar mass, also in solar units, the same as Eqs. in 
what follows. The fitting coefficients arc; seen in Tab. 1, 
the same as Eqs. (5), (6) and (14). And stellar Teg can be 
obtained from L = A-KB^aT^f^. Hence, we can achieve the 
relationship between stellar mass and T^a at ZAMS, seen 
in Fig. 1. 

3.1.2. Stellar luminosity and radius at TMS 

For stellar luminosity and radius at TMS, we modify the 
former shapes (Hurley et al., 2000) and fit again, and the 
expressions are 



anM^-^ + aisM^ + aigM* 



020 



7Vfa23+1.8 



021 + a22M5 + M«23 

rlO+rllM''^^+rl3M''^'' 

-Rtms = { 



M < 1.4, 



M > 1.4, 



(5) 



(6) 



12 3 4 

Fig. 1. The relationship between stellar mass and Tes at 
ZAMS. Points are stellar data calculated by Eggleton's 
code, and the line is taken from fitting formula, the same 
as Fig. 2. 



where both Xtms and Rtms are in solar units. 
3.2. HZs of host stars 

The inner HZ boundary is determined by the loss of wa- 
ter via photolysis and hydrogen escape. And the outer HZ 
boundary is determined by the condensation of CO2 crys- 
tals out of the atmosphere (von Bloh et al., 2007). Jones 
et al. (2006) gave the HZ fiux at both inner and outer 
boundaries, as a function of T^g. 



Sin 
Sout 



(4) ^0 



4.190 X IQ-'^T^g - 2.139 x IQ-^Tefi + 1.296, (7) 



= 6.190 X 10-"T4f - 1.319 X lO'^Teff -|- 0.2341, (8) 



where Sq is solar constant and Tcs is in Kelvin. As we 
have obtained the relationship between stellar mass and 
Teff, we can obtain the HZ boundary fiux of host stars 
with masses from 0.08 to 4.00 Mq, seen in Fig. 2. The 
more massive is the host star, the higher is T^g, and the 
more is the HZ fiux at both inner and outer boundaries. 

As L = AttcPS, and the three parameters arc just as 
Lq, AU and ^0, for our Earth. Hence, the distances at 
both inner and outer HZ boundaries are given by 



din 

AU 

dout 
'AU 



Sin/ Sq 

,Sout/SQ 



-,1/2 



1/2 



(9) 



(10) 



According to the data of Fig. 2 and Eqs. (3), (9) and (10), 
we can obtain the distances of HZs around host stars with 
masses from 0.08 to 4.00 Mq at ZAMS, seen in Fig. 3. 
Using the same method, we achieve the distances of HZs 
around host stars at TMS, also seen in Fig. 3. 
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M [Afg] 

Fig. 2. The HZ flux at both inner boundary (upper line) 
and outer boundary (lower line), considering the Tcs ef- 
fect. 




Fig. 3. The HZs of host stars at ZAMS and at TMS, taken 
the Teff effect into consideration. The lines with points 
denote the HZ boundary distances at ZAMS, and the lines 
with asterisks denote the HZ boundary distances at TMS. 
Points and asterisks are stellar data taken from Eggleton's 
code, and the lines are taken from fitting formulae. 



3.3. Probability distribution of terrestrial planets in 
HZs around host stars 

We take that terrestrial planet is the planet with mass 
from 0.3 to 10 Earth masses (Williams et al., 1997; Menou 
& Tabachnik, 2003; Valencia et al., 2006), and semimajor 
axis is the distance between the planet and its host star. 
Using the data of Fig. 2a in Ida & Lin (2005), we count 
the number of terrestrial planets whose semimajor axis lie 
in the HZs of the host stars at ZAMS, noted as np{M). 
The masses of host stars are 0.20, 0.40, 0.60, 1.00 and 
1.50 Mq, respectively. And we denote the probability of 
terrestrial planets in HZs around host stars as p{M). As 
the evolution of a total of 20 000 planets is calculated for 
each run, the relationship between p{M) and np{M) is 



Table 2. The distances at both inner and outer HZ 
boundaries, and the values of np{M) and p(M), for host 
stars with masses 0.20, 0.40, 0.60, 1.00 and 1.50 Mq 



Stars 


din 




np{M) 


p{M) 


0.20 


0.081288 


0.162912 


52 


0.064707 


0.40 


0.148050 


0.295742 


87 


0.108259 


0.60 


0.248416 


0.494705 


312 


0.388241 


1.00 


0.697398 


1.391183 


745 


0.927049 


1.50 


1.525654 


3.085378 


794 


0.988023 



In order to obtaining the value of kp, we need the 
data about planets in our solar systems. The number of 
planets with masses more than 0.1 Earth mass, for the host 
star with mass 1.00 Mq, is 6429, according to the data of 
Fig. 2a in Ida & Lin (2005). And there are 8 planets with 
masses more than 0.1 Earth mass, in our solar system. 
Hence, we can get the value of kp from 



= kp X 6429/20000, 



(12) 



so that kp = 24.88723. Therefore, we get the linear rela- 
tionship between p{M) and np{M), just as 



p{M) = 24.88723 x np(M)/20000 

= 1.2443615 X 10"^ x np{M). 



(13) 



For the host star with mass 0.20 Mq, the distances 
at both inner and outer HZ boundaries are 0.081288 and 
0.162912 AU, respectively. And the np{M) is 52, so that 
p(M) = 0.064707, according to Eq. (13). Using the same 
method, we also get the np{M) and the p(M) for host 
stars with masses 0.40, 0.60, 1.00 and 1.50 Mq, seen in 
Tab. 2. According to the p{M) for these five stars, we 
obtain p(M) around host stars with masses from 0.08 to 
4.00 Mq, seen in Eq. (14) and Fig. 4. 



p(M) = { 



ogeM'^^^ + a38, 
a39(M - a4o)°-^ + 041, 



M < 0.60, 
M > 0.60. 



(14) 



It may be too rough that p{M) is obtained by extrap- 
olation, with stellar mass more than 1.50 Mq. However, 
it is seen that p{M) increases with increasing stellar mass 
and the increasing trend becomes slower and slower, for 
host stars with masses from 0.60 to 1.50 Mq. And p{M) 
lasts this increasing trend, for host stars with masses more 
than 1.50 Mq. In addition, the value of ^(M) for stars with 
masses more than 1.50 Mq in the Milky Way is very low. 
Therefore, it is acceptable that the formula is fitted from 
five points, although it is not very perfect. 

3.4. Probability distribution of terrestrial planets in 
HZs around host stars in the Milky Way 

We use Prob{M) to denote the probability of terres- 
trial planets in HZs around host stars in the Milky Way. 
Therefore, Prob{M) is equal to p{M) multiplied by ^(M) 
seen in Eq. (1). 



p{M) = kp X ?7p(M)/20000. 



(11) Prob{M) = p{M) X ^(M), 



(15) 
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12 3 4 

Fig. 4. The probability distribution of terrestrial planets 
in HZs around host stars. The data of five points is ob- 
tained from the distribution of semimajor axis and mass 
of planets (Ida & Lin 2005) and our results about HZs of 
host stars at ZAMS. The data of line is fitted from the 
five points. 




Fig. 5. The probability distribution of terrestrial planets 
in HZs around host stars in the Milky Way. 



also intuitively seen in Fig. 5. And Prob{M)dM is the 
probability of terrestrial planets in HZs around all the 
host stars with masses between M and M + dM in the 
Milky Way. 

When stellar mass is less than 0.35 Mq, p{M) increases 
gently, but £,{M) decreases dramatically. Hence, Prob{M) 
decreases with increasing stellar mass. When stellar mass 
is from 0.35 to 0.62 M©, p{M) increasing rate is more than 
^(M) decreasing rate, so that Prob{M) increases with in- 
creasing stellar mass. And Prob{M) meets a peak value, 
with stellar mass 0.62 Mq. When stellar mass is more 
than 0.62 M©, (,{M) decreasing rate is more than p{M) 
increasing rate, so that Prob{M) decreases with increasing 
stellar mass. 



3.5. The number of terrestrial planets in the Milky 
Way 

The number of stars in the Milky Way is about 300 billion. 

And observations tell us that upwards of 50% of the stars 
populating galaxies are expected to be in binary or higher- 
order multiple systems (eg., Duquennoy & Mayor, 1991; 
Richichi et al., 1994). Previously, it was generally thought 
that these systems are not suited to life survive. Recently, 
it was pointed out that terrestrial planets may be stable 
in binary systems for a long time (eg., David et al., 2003; 
Fatuzzo et al., 2006; Saleh & Rasio, 2009). Hence, the 
number of host stars that may be suitable for life survive 
in the Milky Way is about 200 billion, noted N as the 
follows. Therefore, the number of terrestrial planets in HZs 
of host stars in the Milky Way is 



planet 



= N / i 

Jo.08 



p{M)^{M)dM = 4.5507 x 10 



10 



(16) 



As the mass range of M, K, G and F type stars are 
about from 0.08 to 0.50 Mg, from 0.50 to 0.80 Mq, from 
0.80 to 1.05 Mq and from 1.05 to 1.40 Mq (Grie^meier 
et al., 2005), respectively. Therefore, the number of ter- 
restrial planets in HZs for M, K, G and F type stars are 
11.548 billion, 12.930 billion. 7.622 billion and 5.556 bil- 
lion, respectively. 

These mean that the number of terrestrial planets 
around K type stars is the most in the Milky Way. Because 
both the number of K type stars and the p(M) around K 
type stars are great. And the p(M) around M type stars 
is very low, although the number of M type stars is huge. 
On the contrary, the number of G type stars is minor, 
although the p{M) around G type stars is high. 

In order to estimating the number of terrestrial planets 
with lives or even intelligent lives in the Milky Way, we 
use the critical step model of biological evolution (Cater 
& McCrea, 1983; Watson, 2008). The model assumes that 
biological evolution needs sequentially pass several criti- 
cal steps. And the successful probability for each step is 
very low, eg., 10% (Watson, 2008). The five step model is 
perfect, as it matches well with the observed date. Watson 
(2008) pointed out that the five critical steps are life emer- 
gence (primary life), prokaryotes to eukaryotes, asexual 
clones to sexual populations, cell differentiation (complex 
life) and primate societies to human societies (intelligent 
life). 

As the number of terrestrial planets in HZs of host 
stars in the Milky Way is about 45.5 billion, and biological 
evolution needs enough time. Therefore, the number of 
terrestrial planets with primary lives in the Milky Way is 
about 4.3 billion. And the number of terrestrial planets 
with complex lives and intelligent lives in the Milky Way 
are about 3.7 million and 0.36 million, respectively. 

4. Discussion 

There are many planets with semimajor axis about 0.03 
AU around host stars with different masses, most of which 
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are terrestrial planets (Ida & Lin, 2005). And the HZs of 
host stars with masses about 0.10 Mq is approximately 
just as 0.03 AU. Hence, there must be many terrestrial 
planets in the HZs for these stars. Considering that there 
is huge number of stars with masses about 0.10 Mq, the 
probability of terrestrial planets in HZs around these stars 
must be very high in the Milky Way. 

In the Monte Carlo simulations of planet formation 
processes before 2008 (eg., Ida & Lin, 2005), it neglected 
the effects of type I migration of protoplanetary embryos 
due to their tidal interaction with their nascent disks. And 
it is found that the type I migration provides a self-clearing 
mechanism for planetesimals in the terrestrial planet re- 
gion (Ida & Lin, 2008), which may modulate the distribu- 
tion of semimajor axis and mass of planets. 

Ultraviolet radiation is a double-edged sword to life. 
If it is too strong, the terrestrial biological systems will 
be damaged. And if it is too weak, the synthesis of many 
biochemical compounds can not go along (Buccino et al., 
2006). In our future work, we are also plan to study the 
ultraviolet habitable zones of host stars. 

5. Conclusion 

Firstly, we obtain the HZs of host stars with masses form 
0.08 to 4.00 Mq at ZAMS and at TMS, considering the 
Teff effect. Secondly, we give the probability distribution of 
terrestrial planets in HZs of host stars. Thirdly, we calcu- 
late the number of terrestrial planets in HZs of host stars 
is 45.5 billion in the Milky Way, and find that the number 
of terrestrial planets around K type stars is the most in the 
Milky Way. Finally, we present discussions about the host 
stars with masses about 0.10 Mq and the effects of type I 
migration, and introduce our future work. One may also 
send any special request to guojianpol982@hotmail.com 
or guojianpol6@163.com. 
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